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ABSTRACT 

Primordial stars are key to primeval structure formation as the first stellar components 
of primeval galaxies, the sources of cosmic chemical enrichment and likely cosmic 
reionization, and they possibly gave rise to the supermassive black holes residing at 
the centres of galaxies today. While the direct detection of individual Pop III stars will 
likely remain beyond reach for decades to come, we show their supernova remnants may 
soon be detectable in the radio. We calculate radio synchrotron signatures between 
0.5 — 35 GHz from hydrodynamical computations of the supernova remnants of Pop III 
stars in minihaloes. We find that hypernovae yield the brightest systems, with observed 
radio fluxes as high as 1 — 10 /iJy. Less energetic Type II supernovae yield remnants 
about a factor of 30 dimmer and pair-instability supernova remnants are dimmer by 
a factor of more than 10,000. Hypernovae radio remnants should be detectable by 
existing radio facilities like eVLA and eMERLIN while Type II supernova remnants 
will require the Square Kilometre Array. The number counts of hypernova remnants 
at z > 20 with fluxes above 1 /iJy are expected to be one per fifty square degree field, 
increasing to a few per square degree if they form down to z = 10. The detection of 
a z > 20 Type II supernova remnant brighter than 1 nJy would require a 200-300 
square degree field, although only a 1-2 square degree field for those forming down to 
z = 10. Hypernova and Type II supernova remnants are easily distinguishable from 
one another by their light curves, which will enable future surveys to use them to 
constrain the initial mass function of Pop III stars. 

Key words: radiation mechanisms: non-thermal - supernovae: general - galaxies: 
formation - cosmology: theory - radio continuum: galaxies 
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1 INTRODUCTION 

The cosmic Dark Ages ended with the formation of the first 
stars in 10 5 — 10 6 Mq cosmological haloes at z ~ 20 — 30. 
Primordial (or Pop III) stars are the key to understanding 
primeval galaxies, the onset of early cosmological reioniza- 
tion and chemical enrichment, and the origins of the su- 
permassive black holes found in most massive galaxies to- 
day. Unfortuna tely, in spite of their extreme luminosities 
i|Schaererl l2002) and the arrival of next-generation near in- 
frared (NIR) observatories such as the James Webb Space 
Telescope (JWST) and the Thirty-Meter Telescope (TMT), 
individual Pop III stars will remain beyond the reach of di- 
rect detection for decades to come. For now, there are no ob- 
servational constraints on either their masses or their rates 
of formation. 

On the simulation frontier, there has been a 
gradual shift in paradigm over the past decade from 
single 30 — 300 Mq stars forming in isolation in 



Scottish Universities Physics Alliance 



jBromm et all |2002|: 



Nakamura fc Umemura 
2008: 



2007, 



haloe s 

200ll: lAbel et all 12002) ; lO'Shea fc Normanl .. . 

Wise fc Abelll2007t ) to binaries dTurk et al.ll2009F and more 

recently to the possibility of 20 — 40 Mp, stars forming 

in small multiples of up to a dozen (IStacv et alj 
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Clark et al 



Stacy et al 



20111 : iGreif et all 1201 ll ; iHosokawa et all 



20101; 
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20 111 ; I Greif et al.ll2012l ). However, these simula- 



tions do not estimate Pop III stellar masses by modeling the 
formation and evolution of the stars. Instead, they usually 
derive them by comparing infall rates at the center of the 
halo at very early stages of collapse to Kelvin-Helmholz 
contraction times to place upper limits on the final mass of 
the star. No simulation currently bridges the gap between 
the initial formation and fragmentation of a protostellar 
disk and its photoevaporation up to a Myr later, so there 
are no firm numeri cal constraints on the Pop III IMF either 
(see IWhaler]|2012l . for a recent review of primordial star 
formation) . 

Pop III stars profoundly transformed the haloes that 
gave birth to them, dri ving out most of their baryons in su- 
personic ionized flows (| Whalen et al.l [20041 ; iKitavama et alj 
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|2004 lAlvarez et all 120061 ; I Abel et all 120071') and later 
exploding as supernovae (SNe) dBromm et al l 120031; 
Kitavama fc Yoshidall2005l; iGreif et al.ll2007l: IWhalen et al l 



20081 ). Ionization fronts from these stars also engulfed nearby 



haloes, either promoting or suppressing star formation in 
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; Hasegawa et al.l 20091; Susa et al. 


2009 
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). As each halo hosted consecutive 



cycles of stellar birth, H II region formation and SN ex- 
plosions (jYoshida et al.l [2007) , gravity and accretion flows 



merger time-scales of r 



20 Myr at z ~ 20 (e.g. 


Johnson et al.l 


1201(4 IWise et alj|2012 


). At the same 



time, the first SNe enriched the IGM with metals and dust, 
triggering a transition from Pop III to Pop II star for- 
mation (e. g. Macke v et al. Il2003l; ISmith fc Sigurdssonll2007l; 



ISmith et al.ll2009l ; lNagakura et alj|2009l ; IChiaki et al.ll2012h . 

How these two processes determined the masses and forma- 
tion rates of stars in primeval galaxies at 2 ~ 10 - 15, and 
therefore their spectra and luminosities, is unknown. 

Stellar evolution models show that the final fates 
of the first stars rest ed upon their masses at birth 
(|Heger fc Wooslevl 120021): 15 - 40 M B P op III stars died in 
core-collapse SNe jjoggerst et~al] l20ld ) and 40 - 140 M 
stars collapsed to black holes, p erhaps with violent pulsa- 
tional mass loss prior to death |Wooslev et al.l [2007V Pop 
III stars from 140 to 260 Mq died as pair instability (PI) 
SNe, extremely energetic thermonuclear explosions with en- 
ergies of up to 100 those of Type la and Type II SNe (Jog- 
gerst & Whalen 2011; Chatzopoulos & Wheeler 2012 extend 
this lower limit down to 65 Mq for rotating stars). Some 
40 — 60 Mq Pop III stars may have died as hypernovae, 
with energies intermediate to those of core-collapse and PI 
SNe l|lwamoto et al.ll2005l ). Attempts have been made to 
indirectly constrain the masses of the first stars by recon- 
ciling the cumulative nucleosynthetic yield of their super- 
novae to the chemical abundances fou nd in ancient, dim 
meta l -poor stars in the Galactic halo (|Beers fc Christliebl 
120051 ; iFrebel et all [2005). some of which may be contami- 
nated with the ashes of the first generation. For example, 
Ijoggerst et all l|2010l ) recently discovered that the average 
yields of 15 — 40 Mq Pop III SNe agree well with the fos- 
sil abundance s measured in ~ 130 stars with m etallicities 
Z < 10~ 4 Z @ l|Cavrel et al J 12004 lLai et al.l 120081 1 . suggest- 
ing that low-mass primordial stars may have synthesized 
most metals at high redshift. However, stellar archaeology 
is still in its infancy because of small sample sizes, systemat- 
ics in the measurements of some elements, and because the 
imprint of metals from first-generation stars on the second 
is not well understood. 

The direct detection of Pop III SNe may be the best 
prospect for probing the earliest generation of stars in the 
near term. Primordial supernovae may be 100,000 times 
brighter than their progenitors, or, at slightly lower red- 
shifts, the primitive galaxies in which they reside. Their 
transience easily distinguishes them from early galaxies, 
with which they otherwise overlap in colour-colour space. 
Previous studies have investigated detection limits for PI 

1 ; " — 11 1 J 1| -r 

SNe at z ~ 6 ( Scannapicco et al. 2005: Tanaka et al. 2012 , 
for 6 < z < 15 jPan et al.ll201ll ; lMoriva et alj|2012l . Whalen. 
Joggerst et al. in preparation; Whalen, Even et al. in 



preparation), and in very approximate terms for z ~ 30 
l|Hummel et al.l |2012| ). Whalen, Fryer et al. (in prepara- 
tion) and Whalen, Frey et al. (in preparation) show that 
JWST will detect PI SNe beyond z ~ 30 and that the 
Wide Field Infrared Survey Telescope ( WFIRST) and the 
Wide-field Imaging Surveyor for High-Redshift ( WISH) will 
detect them out to z ~ 15 - 20 in all-sky NIR surveys. 
Unfortunately, it may be a decade before such observations 
are possible, given that JWST and WFIRST will not be in 
operation before 2018 and 2021, respectively. 

In the meantime, it may be possible to discover the first 
generation of supernovae in cosmic backgrounds. Pair insta- 
bility SNe deposit up to half of their en ergy into the CMB via 
inverse Compton scattering at z ~ 20 jKitavama fc Yoshidal 
120051 ; IWhalen et afll200Sl ) . lOh et all |2003l ) have found that 
Pop III PI SNe may impose excess power on the CMB on 
small scales via the Sunyayev-Zeldovich effect. Primordial 
SNe may also be manifest as fluctuations in the NIR back- 
ground because they are so much more luminous than their 
host protogalaxies at high redshift. In this paper, we show 
Pop III supernovae may also be revealed through the de- 
tection of the radio emission from their remnants by both 
current and future radio observatories such as the eVLAn, 
eMERLII^, ASKAF0 and the Square Kilometre ArrajQ 
(SKA). 

Previous studies of the observability of collapsing struc- 
tures at z > 10 in the radio have centered on de- 
tecting the 21c m signal of minihaloes during the cosmic 
dark ages (e.g. iFurlanetto fc Loebl |2002| ; llliev et all l2002j ; 
IShapiro et al.l l2006; Meiksin 2011), not cosmic explosions at 
first light. There are several advantages to searching for the 
first supernovae through their radio synchrotron emission. 
Firstly, if only certain types of explosions contribute a radio 
signal, this may be used to place limits on the masses of 
their progenitors. Secondly, supernova (and hence star for- 
mation) rates may possibly be determined from the strength 
of the radio signals and the number counts of the sources. 
Finally, depending on the properties of this signal, it may 
be possible to pinpoint the redshift of the explosions. In this 
paper we calculate the radio signal produced by the rem- 
nants of core-collapse SNe, hypernovae and PI SNe to assess 
the detectability of the remnants by current and future ob- 
servatories. We also determine if this signal may constrain 
the masses and formation rates of Pop III stars in high red- 
shift protogalaxies. In §[2] we review how Pop III SN explo- 
sions generate synchrotron emission that could be directly 
detected by current and planned radio telescopes. In §[3] we 
describe the expected radio light curves and number counts 
of the Pop III SN radio remnants. We summarise our conclu- 
sions in §21 Formulae are provided for synchrotron emission 
and self-absorption in terms of the energy densities of the 
relativistic electrons and the magnetic field in an appendix. 

For numerical cosmological estimates, we adopt fl m = 
0.27, fl v = 0.73, H = lOOkms^Mpc" 1 , cr 8h _i = 0.81 
and n = 0.96 for the total mass and vacuum energy 
density parameters, the Hubble constant, the linear den- 
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sity fluctuation amplitude on a scale of 8/i~ 1 Mpc and 
the spectral index, respectively, consistent with the best- 
fitting cosmological parameters for a flat universe as con- 
strained by the 5-yr Wilkinson Microwave Anisotropy Probe 
( WMA P) measurements of the Cosmic Microwave Back- 
ground (|Komatsu et alj|2009h . 



2 THE SYNCHROTRON SIGNATURE OF 
SUPERNOVAE IN MINI-HALOES 

We begin by estimating the expected synchrotron flux from 
a supernova remnant. For relativistic electrons in a mag- 
netic field of energy density u b , the bolometric emissivity of 
synchrotron radiation is 



Ebol ~ UBT Bt 



(1) 



where r se is the characteristic energy scattering time-scale 
for electrons of energy density u e , defined by 

r sc = ™=^-. (2) 
ca T u e 

Here <tt is the Thomson cross section for electron scattering 
and m e is the mass of an electron. 

Allowing for a fraction f e of the thermal energy to go 
into relativistic electrons, so that u e = /eiith for a ther- 
mal ener gy density u t h, and further assuming equipartition 
ub — u„ (|Chevalierlll982r ). the bolometric synchrotron power 
emitted by a single halo is 



PZT * 



car 



drr 2 / e (r) 2 u t 2 h F(p e ;7i,7 u ), (3) 



where the function F(p e ; Ji,fu), described in the Appendix, 
allows for a power-law energy distribution of index p e for 
the relativistic electrons, and 7; and y u are the lower and 
upper ranges of the relativistic 7 factors for the electrons. We 
generally assume conservatively 71 = 1, as non-relativistic 
electrons will not contribute much synchrotron emission. 

We illustrate the expected flux using the computa- 
tion for a 40 Mp hypernova in a 1.2 x 10 7 Mq halo at 
z = 17.3 (|Whalen et al.ll2008l V The synchrotron power may 
be estimated from the post-shock thermal energy density 
of the supernova remnant, which exceeds a few erg cm -3 
over a region around a tenth of a pars ec across durin g 
the first ~ 10 yrs. For f e = 0.01 (e.g. IChevalierj fl982l) . 
this corresponds to a bolometric synchrotron power out- 
put of PZi C ^ 5 x 10 42 ergs- 1 for 7 „ = 300 and p = 2. 
The spectrum will extend to a lower characteristic cutoff 
wavelength, assuming an isotropic magnetic field, of A c = 



(2 7 r 1 '' 2 /3)(m e c 2 /e)7; 



-1/2 



2014 7 ; 



-1/2 



cm « 0.2 cm, 



or upper cutoff frequency u c ~ 130 GHz, corresponding to a 
characteristic specific power at v c of 2 x 10 31 ergs -1 Hz -1 . 
For a source at z = 17.3, this corresponds to an observed 
radio flux of ~ 10 /uJy. The flux will vary with frequency as 
„-(pe-i)/2 _ j,-V2_ " 

The radio spectrum computed 4.7 yrs after the explo- 
sion is shown in Fig. \T\ The free-free radiation is shown as 
well, but is generally much smaller than the synchrotron. 
The energy density in relativistic electrons is sufficiently 
high that synchrotron self-absorption (SSA) is significant at 
the lower frequencies. Allowing for SSA severely attenuates 
the spectrum at v £ 20 GHz, as shown in Fig. [1] 




10 s 10" 
v (Hz) 



Figure 1. Radio power from a 40 Mq hypernova in a 1.2xlO 7 M 
minihalo at z = 17.3 (for p e = 2, -y u = 300 and f e = 0.01). Shown 
are total (synchrotron and free-free) (solid lines) and free-free 
(short-dashed lines) powers. The upper curve (magenta) corre- 
sponds to a model with no attenuation. The middle curves allow 
for synchrotron self-absorption (blue) and free-free absorption as 
well (cyan). The lowest curve (black) further adds hypothesized 
plasma synchrotron attenuation effects. 



Adding free-free absorption degrades the spectrum fur- 
ther at low frequencies. Free- free absorption was neglected 
from the fluid zone at the shock front, where the gas is be- 
coming ionized but has not yet reached the post-shock tem- 
perature. Since the shock front is unresolved, the substantial 
free-free absorption from it is likely greatly over-estimated. 
For a shock front width on the order of the Coulomb mean 
free path, the free-free absorption from the shock front be- 
comes negligible. It will therefore generally not be included 
here. It is remarked, however, that if cool circumstellar gas 
mixes in with the post-shock ionized gas, free-free absorp- 
tion could be non-negligible. Such an effect is beyond the 
capacity of a spherically symmetric code to reproduce, and 
would be difficult to resolve in any case. 

The Razin-Tsytovich plasma effect on the synchrotron 
radiation mechanism may fur ther degrade the escap- 
ing synchrotron rad i ation (e.g. IScheuer fe Williams! Il968l ; 
iRvbicki fe Lightmanl fl98r? ) . We approximate the effect by 
cutting off the production of synchrotron radiation by the 
factor exp( — Vwr/v), where 



vrt = luVp ( 1 + 7«^p ) . 



(4) 



and Up is the electron plasma frequency. This strongly cuts 
off most of the radiation below ~ 0.1^ c , as shown in Fig. [1] 
As the occurrence of the hypothesized effect is unclear, we 
conservatively neglect it further here, but note that the de- 
tection of such greatly diminished low-frequency emission 
would support its reality. 

Radio observations of supernova remnants show a range 
in values of p e , typically 1.5 < p e < 3, and do not strongly 
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Figure 2. Radio power P v (erg s — 1 Hz -1 ) as function of emission frequency u (Hz) from a 40 Mq hypernova in a 1.2 X 10 7 Mg minihalo 
at z = 17.3 (for f e = 0.01). The relativistic electron energy power-law index varies over 1.5 < p e < 3.0 across the panels. (Note the 
change in scale for the p e = 3.0 panel.) The power is shown for 7„ = 50 (dotted lines), 100 (solid lines), 300 (short-dashed lines) and 
1000 (long-dashed lines). 



constrain y u and / e , althoug h values of j u > 100 and 
fe ~ 0.001 — 0.2 are indicated (|Weiler et al. 1986; IChevaiier1 
ll99Sl : ISoderberg et al.|[200r}| . The parameters p e and y u are 
allowed to vary in Fig. [2]for f e = 0.01. The power at v — v c 
is proportional to j't ■ For a given 7„, the peak power de- 
creases with increasing p e . The peak power for a given p e 
varies little with 7„ for "/ u > 100, while increasing "/ u further 
increases the synchrotron cut-off frequency and so extends 
the power-law spectrum to higher frequencies. For p e ~ 3, 
a substantially larger value of either f e or the minimum 7 
factor for the relativistic electrons is required to reproduce 
the same power as for f e < 2.5. 



It is noted that for f e — 0.01, the total synchrotron en- 
ergy radiated will match the energy in relativistic electrons 
after a few years. Either 7 U must then decrease with time or 
f e must be lower or decrease with time (most likely associ- 
ated with a steepening of the power-law energy distribution 
of the relativistic electrons), unless electrons are continu- 
ously accelerated to relativistic velocities by the shocks, as 
is believed to occur in supernova remnants in the Galaxy. 
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Figure 3. Radio light curves for a 40 Mq hypernova in a 1.2 X 
10 7 Mq minihalo at z = 17.3 (for p e = 2, ■ju = 300 and f e = 
0.01). The observed fluxes are shown for bands: 0.5 (dotted), 1.4 
(solid), 3 (short-dashed), 10 (long-dashed), 25 (dot short-dashed) 
and 35 (dot long-dashed) GHz. 



Figure 4. Radio light curves for a 15 Mq Type II supernova in 
a 1.2 X 10 7 Mq minihalo at z = 17.3 (for p e = 2, -y u = 300 
and f e = 0.01). The observed fluxes are shown for bands: 0.5 
(dotted), 1.4 (solid), 3 (short-dashed), 10 (long-dashed), 25 (dot 
short-dashed) and 35 (dot long-dashed) GHz. 



3 OBSERVATIONAL PREDICTIONS 



Deep radio surveys are currently able to reach the few fxJy 
level. The Very Large Array has achieved an rms noise level 
of 1.5 uJ v at 8.4 GHz (jFomalont et aT1l2002l ) and 2.7/iJy at 
1.4 GHz (jOwen fe Morrisonll2008l). Similar levels have been 
reached using eMERLIN (jMuxlow et al.ll2005l ). The Square 
Kilometre Array is projected to be a factor 100-1000 more 
sensitive. 

In Figs [3] and [4l we show the light curves for a 40 M 
hypernova and 15 Mq Type II supernova, respectively, in a 
1.2 x 10 7 Mq minihalo at z — 17.3 in selected bands planned 
for the SKA. The bands lie within the eVLA frequency 
range. The eMERLIN L-band includes 1.4 GHz, while the 
C-band flux light curve would lie between those for the 3 
and 10 GHz bands. Time-dilation produces a gradual rise 
extended over a 100 yr time-scale with an even slower de- 
cline. For the 40 Mq hypernova, a flux exceeding 1 /xjy is 
expected to persist over a span of 300 yrs, with a peak flux 
of ~ 10 /iJy. Most notable is the relatively late rise of the 
flux in the 500 MHz band; it eventually overtakes the higher 
frequency fluxes for a period of 200 yrs while above 1 /uJy. 
A similar trend is found for the 15 Mq supernova, although 
over the much shorter time-scale of decades and with flux 
values lower by 1.5 dex. We also computed the emission 
following a 260 M pair instability supernova, but found it 
so efficiently sweeps through the minihalo that the resulting 
observed radio flux is at most ~ 0.1 njy. Similarly, we found 
the observed radio fluxes from the remnants produced in a 
2 x 10 6 Mq halo by all three supernova progenitor masses 
fell well below 1 njy. 

Allowing for the Razin-Tsytovich effect much reduces 
the flux in the 500 MHz band, as shown in Fig. [5] While it 
continues to rise with time, it lies about 1 dex below the flux 



40 M hypernova 
Razin-Tsytovich effect included 




10 100 
time (yrs) 



Figure 5. Radio light curves for a 40 M hypernova in a 1.2 X 
10 7 Mq minihalo at z = 17.3 (for p e = 2, 7 U = 300 and f e = 
0.01). The Razin-Tsytovich effect is included. The observed fluxes 
are shown for bands: 0.5 (dotted), 1.4 (solid), 3 (short-dashed), 
10 (long-dashed), 25 (dot short-dashed) and 35 (dot long-dashed) 
GHz. 
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at 1.4 GHz for 350 yrs. When it does overtake the 1.4 GHz 
flux, its value lies well below 1 fj, Jy. Comparison with Fig. [3] 
suggests the measured flux in the 500 MHz band may be 
used to detect the presence of the Razin-Tsytovich effect. 

The number of sources visible per year is substan- 
tial. We estimate this fro m the halo collaps e rate using 
the halo fitting function of iReed et al. (2007), adapted to 
the 5 yr WMAP cosmolog ical parameters for a flat uni- 
verse l|Komatsu et al.ll2009T ). The observed formation rate 
of haloes within a solid angle 8£l with a comoving number 
density n^ om (z) for halo masses exceeding Mh, integrated 
over zi < z < Z2, is then 



Nt s (> M h ) 



dz »ifi om (l + z) 



dz 1 



22 rln COIn 

dz ^h—daarfSSl 
dz 



+2 



7 com 

dz rih c-. 



a r 



n^° m (2i)c[aor(zi)] 2 (5ri, 



where aor(z) = (c/Ho) J Q Z dz [f2 m (l + 2) 3 + fi v ] _1//2 is the an- 
gular size distance for a flat universe with present day mass 
and vacuum energy density parameters Q m and Q v , respec- 
tively, and Hubble constant Ho- For z\ > 3, (Ho/c)aor(zi) — 
1. 5 + 1. 9[1 — 2/(l + 2i ) ] to better than 2 percent accuracy. 
The approximation h^™ 1 — (dn^° m / dz){dz / dt) was made, 
although this could be modified by allowing for merger his- 
tories. In the final line, the integral in the line previous was 
neglected as was the halo density at z = 22. A factor 1/(1+2) 
has been included to account for time-dilation. 

We find formation rates for haloes M > 10 7 h" 1 Mq and 
2 > 20 of 0.0016 deg" 2 yr-\ and 0.20 deg~ 2 yr^ 1 at 2 > 10. 
These amount to about 70 to 8000 collapsing haloes per year 
over the sky. The rates are comparable to recent predictions 
for the total production rate of Pop III pair instability su- 
pernovae, with progenitor masses bet ween 140 — 260 Mm . 
based on cosmologica l simulations (e.g. iHummel et"aT, I l2012l ; 
I Johnson et"afl l2012h . assuming one PI SN per minihalo 
formed. On the other hand, small scale simulations suggest 
fragmentation may prevent the formati on of Pop HI star s 
sufficiently massive to form a PI SN (e.g. IStacv et aljfeoid ). 
although subsequent mergers of the protostars into more 
massive ones cannot be ruled out. The star formation may 
also be stochastic given the low number of massive stars 
formed, with hypernovae forming before a more massive PI 
SN progenitor is created. 

The effect of multiple supernovae on the gas in small 
haloes is unknown. The heat input of a supernova is suf- 
ficient to unbind the IGM of a minihalo with mass below 
~ 1O 7 M0, while, once the gas cools , it w il l fall back in a 
more massive system l|Whalen et al.l [§008 ; Mei ksinl l201ll ) . 
The remaining massive stars will further replenish the IGM 
through winds, so a gaseous environment may continue to 
persist on which further supernovae may impact, producing 
synchrotron emitting remnants. If each halo gives rise to at 
least one supernova with progenitor mass exceeding 40 Mq, 
its remnant should then be detectable by eVLA or eMER- 
LIN. If each produces one with a progenitor mass exceeding 
15 M0, SKA should be able to detect the remnant. 



An estimate of the expected number counts of radio 
emitting remnants requires modelling the remnants over a 
broad redshift range. In the absence of simulations covering 
a wide range, we use those discussed above, presuming the 
behaviour of the gas is not very sensitive to redshift. This is a 
coarse approximation, as the gas content of the haloes is ex- 
pected to be sensitive to the redshift of their formation. On 
the other hand, since much of the gas into which the super- 
nova explodes is produced by stellar winds, the environment 
of the progenitor may not be very redshift dependent. An 
improved estimate would require a better understanding of 
star formation within primordial haloes. 

The number of sources with observed flux exceeding 
f v is given by modifying equation ((5| to account for the 
duration 8t e ^ (2) the remnant is visible with a flux exceeding 
U: 

rz2 j com 

N(> fu) = - dz ^— c {a rf5tf{z)5Q,. (6) 

The resulting number counts for haloes collapsing at 2 > 20 
shown in Fig. [6] based on one 40 Mq hypernova per col- 
lapsing halo with a mass exceeding 10 7 h~ 1 Mq. The counts 
(^re fairly flat up to 5/iJy. A source brighter than 1 pjy 
should be visible every 50 square degrees, while one source 
brighter than 1 njy should be visible per 10 square degrees. 
At the earliest epochs of star formation, although the rem- 
nants would still have observed peak fluxes brighter than 

1 /iJy, the frequency of their occurrence becomes very small. 
Only one remnant at 2 > 30 brighter than 1 njy is expected 
over the entire sky . 

I Johnson et"aD (|2012h suggest Pop III star formation 
may persist to 2 < 10. This substantially boosts the counts 
to a few per square degree brighter than 1 /iJy, and more 
than 10 per square degree brighter than 1 njy. The sources 
may be identified by their radio colours. Even more defini- 
tive, however, would be a survey campaign with repeat 
multiband observations extended over several years to trace 
the light curves. 

The corresponding number counts based on one 15 Mq 
Type II supernova per collapsing halo are shown in Fig. [7] 
While the fluxes are considerably weaker than for the more 
energetic hypernovae, the numbers are lower as well as a 
result of their shorter durations. To detect the remnants at 

2 > 20, 200-300 square degree fields would be required to 
detect remnants as weak as 1 njy. The counts for 2 > 10 are 
considerably higher, requiring now 1-2 square degree fields 
to detect the supernova remnants. 

We have also examined the radio absorption signature 
of the remnants against a bright background radio source. 
While synchrotron absorption will produce a strong absorp- 
tion feature at rest frame frequencies below 20 GHz, it lasts 
less than ~ 100 yr. The number of detectable features along 
a line of sight will be 



dN{> T„) 
dz 



com / com\2/i . \ r_tCif / \ 



(7) 



where Vq ™ 1 is the comoving radius of the remnant out to 
which the line-of-sight absorption optical depth exceeds t„ 
and St^(z) is the effective duration of the feature. For a 
characteristic comoving radius of ~ 1 pc, dN/dz ~ 3 x 10 -14 
at 2 = 10, so that the feature would be undiscoverable. At 
late times, at t ~ 0.8 Myr, the outflowing gas cools suf- 
ficiently to produce a 21cm absorption feature that would 
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Figure 6. Number counts of radio remnants above a given flux, 
based on radio light curves for a 40 Mq hypernova in 1.2 X 10 7 Mq 
minihalocs (for p e = 2, 7„ = 300 and f e = 0.01) at z > 20 (upper 
panel) and z > 10 (lower panel). The observed counts arc shown 
for bands: 0.5 (dotted), 1.4 (solid), 3 (short-dashed), 10 (long- 
dashed), 25 (dot short-dashed) and 35 (dot long-dashed) GHz. 
The counts at 25 and 35 GHz would be close to those at 10 GHz 
for 7„ > 1000. 



Figure 7. Number counts of radio remnants above a given flux, 
based on radio light curves for a 15 Mq Type II supernova in 1.2 X 
10 7 M minihaloes (for p e = 2,y u = 300 and f e = 0.01) at z > 20 
(upper panel) and z > 10 (lower panel). The observed counts are 
shown for bands: 0.5 (dotted), 1.4 (solid), 3 (short-dashed), 10 
(long-dashed), 25 (dot short-dashed) and 35 (dot long-dashed) 
GHz. The counts at 25 and 35 GHz would be close to those at 
10 GHz for y u > 1000. 



be measurabe by SKA. We find a characteristic 21cm ab- 
sorption doublet arising from the cooling shell would form 
along lines of sight within ~ 5 pc (proper) of the centre of 
the minihalo. The feature would survive about 1 Myr before 
transforming into a singlet absorption feature indistinguish- 
able from those expected from minihaloes. This corresponds 
to dN/dz ~ 1CF 6 , again too small to be discovered as there 
would still not be nearly an adequate number of bright back- 
ground radio sources to have a fair chance of seeing even one 
feature. 



4 CONCLUSIONS 

We estimate the radio signatures of Pop III supernovae in 
~ 10 7 Mq minihaloes, sufficiently massive to retain their 
baryons and form supernova remnants, based on hydrody- 
namical computations of 15 M Q Type II supernovae, 40 M 
hypernovae and 260 Mq pair instability supernovae at z = 
17.3 within the haloes. We model the synchrotron emission 
and absorption assuming a power-law distribution of rela- 
tivistic electron energies with a total energy density pro- 
portional to the thermal energy of the ionized gas, and 
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equipartition between the relativistic electron and magnetic 
field energies. Allowing for a relativistic electron component 
with total energy one percent of the gas thermal energy 
is sufficient for a hypernova to produce observable fluxes 
at 0.5 — 35 GHz exceeding 1 /xjy, and for the less massive 
Type II supernova to produce observable fluxes exceeding 
10 nJy at 0.5 - 25 GHz. The PI SN expels much of the 
halo gas, leaving behind a supernova remnant that produces 
fluxes of at most 0.1 nJy. A halo mass exceeding 10 8 ft -1 Mq 
would be sufficient to retain the baryons following a PI 
SN. Although we do not currently have computations of 
PI SN in such massive haloes, we may estimate their for- 



-l 



-6 . 



1 for those forming at z > 20. These 
are about a factor 20-200 smaller than the formation rates 
of minihaloes with masses exceeding 10 7 h~ x Mq, but would 
still produce a detectable yield of PI SNe remnants over the 
sky provided their synchrotron fluxes were sufficiently high. 

Hypernovae produce observed light curves that rise 
gradually over a period of 100 yrs, before declining after 
200-300 yrs. By contrast, the Type II curves have rise times 
of 10-20 yrs, and decline abruptly after 30 yrs. The 500 MHz 
flux lags the higher frequency bands for both supernovae and 
dominates the emission at its peak as the synchrotron spec- 
trum reddens due to the weakening of the supernova shock. 
Light curves computed with and without the hypothesized 
Razin-Tsytovich effect show that the observed 500 MHz flux 
is substantially diminished by the effect, suggesting it may 
provide a useful means of testing for its existence. 

Estimating the number counts of the supernovae from 
the minihalo formation rate at z > 20, we find a supernova 
formation rate of about one per 600 square degrees per year. 
If the supernovae are able to form down to z = 10, the rate 
increases to one per 5 square degrees per year. On average, 
one z > 20 hypernova radio remnant with flux exceeding 
1 fiJy should be detectable at any given time per 50 square 
degree field, and one with a flux above 1 nJy per 10 square 
degree field. If Pop III star formation persists to z < 10, 
then a few hypernova remnants brighter than 1 fijy would 
be visible per square degree, and more than ten per square 
degree brighter than 1 nJy. 

Because of their smaller explosion energies, Type II su- 
pernovae will produce weaker radio remnants with shorter 
durations. Fields of 200-300 square degrees would be re- 
quired to detect remnants forming at z > 20 with fluxes 
above 1 nJy. The numbers are much improved for super- 
novae forming at z > 10, requiring instead 1-2 square degree 
fields to detect their radio remnants. 

Pop III supernova remnants with radio fluxes exceeding 
1 (iJy should be visible in radio surveys using the eVLA or 
eMERLIN. Weaker few nJy sources would be detectable by 
the SKA in large numbers. An observing campaign extended 
over several years would be able to follow the light curves, 
yielding invaluable data on the formation of the first stars in 
the Universe and their impact on their local environments. 

The readily distinguishable differences between the ra- 
dio light curves of hypernova and Type II supernova rem- 
nants provide the possibility of exploiting the radio emission 
to infer the initial mass function of the first stars. The fre- 
quency with which these explosions are detected over the 
sky will also constrain Pop III star formation rates through 
cosmic time along with all-sky NIR surveys by WFIRST and 
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Figure Al. Synchrotron function F(p e ', 7i, T«) f° r relativistic 
electron power-law energy distribution over 7; < 7 < -y u with 
index p e . Shown for p e = 1.5 (solid line), p e = 2.0 (dashed line), 
p e = 2.5 (dot-dashed line) and p e = 3.0 (dotted line). (A value 
7; = 1 is adopted.) 



WISH and deep-field surveys by JWST. Their discovery in 
the radio will permit followup by deep-field observations of 
the host galaxies by JWST and the TMT. Surveys detect- 
ing Pop III radio remnants at 10 < 2 < 15 may also con- 
strain the multiplicity of their occurrence within individual 
primeval galaxies, and even the large-scale distribution of 
the galaxies themselves. The detection of primordial super- 
nova remnants will be one of the most spectacular discover- 
ies in high-redshift radio astronomy in the coming decade. 



APPENDIX A: SYNCHROTRON EMISSION 
AND ABSORPTION 

The synchrotron emissivity for a power-law relativistic elec- 
tron energy distribution with index p e extending between 
the 7 factors 7; < 7 < "/ u may be expressed concisely as 



3-Pe tbol 

2 u c 



\ < pe < 3 



(Al) 



(A2) 



l°g("c/Vmin) (v c J !Pe — 3, 

where u c is the synchrotron upper cutoff frequency 

3 7,^eBsina 
v c = 

4-7T m e c 

for pitch angle a, and tboi is the bolometric emissivity 

eboi = UsT S e F(p e ;ji,y u ). (A3) 

Here, Ms is the energy density of the magnetic field, the 
characteristic electron scattering time r se is given by equa- 
tion ([2]) and the function F(p e ; 7;, 7„), shown in Fig. IA1I 
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X < 



f 2-p e ,, l-(7l/7u) 3 ~ Pe 
3- Pe 7« 1 _ (7i/7ii) 2- Pe 

7«T^fe ;Pe = 2 (A4) 



log(7u/7i) 
_,,„, 'og(7u/7l) 



;p e = 3. 



The inverse attenuation length due to synchrotron self- 
absorption may be expressed as 

ot v = r-Q 1 ( — — j ( — ) (Pe G(p e ;7;,7 u ), (A5) 



where ro is the classical electron radius and 



2-Pe 
l-(7!/7u) a - p ': 
1 

log(7u/7l) 



;p e = 2. 



t) 



(A6) 
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